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ABSTRACT 

A knowledge of the complex dielectric constant of soils is essential in the interpretation 
of microwave airborne radiometer data of the Earth’s surface, Measurements were made at 
37 GHz on various soils from the Phoenix, Ariz., area. Extensive data have been obtained on 
e * e f - /e" for dry soil and soil with water content in the range from 0.6 to 35 percent by 
dry weight. Measurements were made in a two-arm microwave bridge and results were cor- 
rected for reflections at the sample interfaces by solution of the parallel dielectric plate 
problem. The maximum dielectric constants are about a factor of 3 lower than those re- 
ported for similar soils at X-band frequencies, 
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DIELECTRIC CONSTANTS OF SOILS AT MICROWAVE FREQUENCIES 


F, E. Geiger and Donald Williams 
Goddard Space Flight Center 


INTRODUCTION 


Purpose of Measurement 

In the microwave region the emergent thermal radiation from a semi-infinite medium is 
proportional to the brightness temperature (by the Ray leigh-Jeans approximation, refs. 1 
and 2) 


T b = (i - |/? 12 | 2 >/ a(z)T(z) exp 

Jo 



a(£) d( 


dz 


( 1 ) 


where |/C 2 | 2 is the free-space power reflection coefficient of the semi-infinite medium 2 with 
respect to the medium 1 at the interface z = 0; 1 - |R 12 | 2 * s ^ le em issivity of medium 2, and 
is related to the thermodynamic temperature T(z = 0) by 1 - |i? 12 | 2 - T b /T( 0); and a(z) is 
the free-space attenuation coefficient of the electromagnetic wave in medium 2 (not to be 
confused with the guide absorption constant a. defined later); |R J2 | 2 is related to the dielec- 
tric constant e 2 of medium 2 by the Fresnel equations. Consequently, the emergent radia- 
tion becomes a direct function of the dielectric constant e 2 . But equation (1) also shows 
that T b is a function of characteristics of the interior of medium 2. Radiometric measure- 
ments of thermal radiation are therefore dependent on a detailed knowledge of surface 
emissivities (1 - R(e) 2 and absorption constants ot(z) and their dependence on physical con- 
ditions in the interior of the medium. 


This paper reports on the determination of the complex dielectric constant e = e' - je” 
of a number of soil samples from the Phoenix, Ariz., area as a function of moisture content 
at a frequency of 37 GHz. 


Review of Existing Methods 

The experimental method for the measurement of e must be suitable to measure the 
real part of the dielectric constant 1 in the range 2.5 < e f /e Q <10 and the imaginary part 
from 0.05 < e"/e 0 < 1 0. 

*SI units are used throughout; e 0 = 8.85 X 10“ 12 F iri' 1 . 
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Four previously used methods were considered for these measurements; 

( 1 ) Truttemishion measurements between two points in an infinite line 

(2) Input impedance measurements on a finite line 

(3) Cavity resonance methods 

(4) Free-space transmission measurements 

A thorough discussion of the First two is given by Westphal (ref. 3), and specific examples of 
application of the first to water are given by Buchanan (ref. 4), to soils and water by Straiton 
and Tolbert (ref. 5), and to ice and snow by Cumming (ref. 6). The second method has 
been used by von Hippel (ref. 7) on a wide variety of materials, solid and liquid, and for a 
wide range of e f and e". 

Cylindrical cavity resolutions have been uud in the transverse magnetic (TM 010 ) and 
transverse electric (TE 01rt ) modes with disk- and rod-shaped samples. Solutions of Maxwell’s 
equations for these various configurations and modes have been worked out by Horner et 
al. (ref. 8), and Tkach (ref. 9). Tkach measured a large variety of solids and liquids with 
relatively high e ' and e" (e'/e 0 ~ 20, e n /e 0 ~~ 2). 

Straiton and Tolbert (ref. 5) have used fiee-space transmission methods on both soil 
and water. 

Choice of Methods 

The suitability of the standing-wave detector for measuring a wide range of impedances 
and the elegance of von Hippefs method (ref. 7) makes the input impedance method very 
attractive. Serious drawbacks for accurate measurements would seem to be the apparent 
unavailability of commercial precision voltage standing-wave ratio (VSWR) bridges in the 
centimeter and millimeter ranges. Von Ilippel’s apparatus for a wavelength of X = 6 cm had 
a probe position accuracy of 0.0025 cm (0.001 in.). Greatest sensitivity is achieved wh°n 
the length of lossy sample is approximately an odd multiple of guide wavelengths (in the 
dielectric). For free-air wavelengths of 0.8 and 1 ,6 cm, optimum sample lengths are found 
to be on the order of 5 and 1 2 mm, respectively. For soil samples this probably would give 
unrepresentative results. Very lossy samples would be insensitive to the length of the sam- 
ple because the sample impedance approaches the guide wave impedance. 

In cylindrical cavity resonance methods at high frequencies, one has to contend with 
small cavities. TM modes have high axial electric fields and result in smaller cavities than 
do TE modes. The latter have no axial fields and may be suitable for relatively high loss 
measurements. Using Tkach’s (ref. 9) equations for TE 01/J modes and coaxial cylindrical 
samples, calculations show that for 19 and 37 GHz, and diameter to length of c dty ratios 
of ( D/L ) 2 ~ 1/10, sample diameters are less than 0.5 and 0.25 cm, respectively. This cor- 
responds to an upper limit in the loss tangent of approximately tan 5 < 0.005 for both fre- 
quencies. These figures show dearly that soil sample measurements are not practical in this 
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frequency range by cavity resonance methods. Considerations for 5 GHz show that speci- 
mens 1 cm in diameter could be used for tan 5 < 0.01. 

Free-space transmission measurements have the advantage that sample preparation is 
very simple for high-loss samples. The sample can be placed on a plastic tray between two 
microwave horns of a bridge iref. 10). For low-loss samples this may be impractical. Fur- 
ther, some factors in the measurements are difficult to assess, such as spherical wave fronts, 
horn antenna patterns with samples of different refractive index, and 1/r effects in large 
samples. 

An elegant free-space ellipsometer method apparently was introduced by Pozdniak 
(ref. 1 1). Elliptically polarized radiation reflected from a surface can be analyzed in terms 
of the polarization coefficient. The polarization coefficient in turn is a function of the 
dielectric constant of the surface. However, radiation incident on a rough surface suffers 
depolarization changing the polarization coefficient (phase and amplitude). Hence this 
method probably has only limited applicability for soils. 

The last method to be discussed is the transmission method on an infinite line, as sug- 
. gested by T. Wilheit (personal communication). Although it is a straightforward measure- 
ment of attenuation and phase, its accuracy is limited primarily by the calibration accuracy 
for attenuation and phase shift. It can, however, cover a wide range of e' and e" with good 
accuracy, especially when used in a bridge configuration. It has the advantage that sample 
sizes are of convenient and practical proportions for soil measurements. As a result, this 
method was chosen for measurements at 37 GHz. 

PRINCIPLES OF A MICROWAVE BRIDGE 

Fundamental Equations 

The attenuation constant a and phase constant j3 for a TE wave in a rectangular guide 
can be written in the form (ref. 1 2) 


where 


7 — a + ffi 


(« + //?) 


2 - 


2tt 


^0 L 



- c 2 /i 0 e'-fn 0 c 2 e" 


e-e* - je” 


( 2 ) 

(3) 

(4) 


7 = propagation constant 
a = attenuation constant, Np irf 1 
f 3 = phase constant, rad nr 1 

X 0 = free space wavelength, m « c/v Q , c = 3 X 10 8 m s~ 1 


N 
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X (1 = aJ2 - cutoff wavelength of guide for TF 1U mode, m, where a is width of guide, m 
= permeability of free space, 1.257 X 50 6 II m' 1 
c' = real part of dielectric constant, F m 1 
c" ~ imaginary part of dielectric constant, F m 1 
e = permittivity of medium, F nr 1 
Other definitions to be used are 






(5) 

( 6 ) 


where 


fi 0 == propagation constant in guide for free space 
j3 c = propagation constant in guide with medium of permittivity e 
\o = wavelength in free space, m 
X^ c = guide wavelength in medium with permittivity e, m 


and where 


\o ~ \) 


\e = \ 



- 1/2 


- 1/2 


Solving equation (3) for e' and e" gives 



(7) 

( 8 ) 


(9) 

( 10 ) 


Method of Determining a and (3 e 

The definitions of ex and /S, suggest their method of measurement. The TE wave ampli- 
tude attenuation in a sample of length L is given by e '~ ocL . If the net power attenuation 
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(attenuation witii sample in g,.ide minus attenuation in the empty guide) is A dB, a can be 
calculated from 

A Np m* 1 (li) 




8 . 6 /, 


Similarly, if the length of the sample is changed by A/, and the accompanying phase 
shift in the transmitted wave is A0, then 

A0 


0 . 


A / 


+ fir 


( 12 ) 


Very lossy samples (losses greater than 30 dB cm"* 1 require different measurement and cal- 
culation procedures for and will therefore be discussed as a separate case. The length oi the 
sample cell for these samples is very small, typically on the order of 0,3 to 0.5 cm. Longer sam- 
ple lengths lead to very v/eak null signals. The short length makes it impractical as well as inaccu- 
rate to determine P e by changing the sample length as outlined previously. As will be shown, the 
accurate determination of/3, depends on the condition A/ «.L ; but this condition now is no 
longer met because Al~L, Consequently, the procedure is changed, and one determines the 
phase shift in the sample of length/,. This shift A 0 is of manageable magnitude (i.e„ 2i r to 3ir 
rad) and can be measured unambiguously. But this A 0 now determines the phase of the trans- 
mitted wave. (See detailed discussion in the section entitled “Reflection Errors.”) However, if 
the amplitude and phase of the wave emerging from the sample are known, /3 tf can be calculated. 

The measurement of A and A0 can be done in a bridge setup in which the wave from 
the sample is nulled by a wave from the same source appropriately attenuated and shifted 
in phase. 

The experimental arrangement to do this can take numerous forms. Buchanan (ref. 4) 
used a hybrid T as the bridge element with attenuator and phase shifter in E and H arm, 
respectively. The sample cell can be put in either arm. The method used here is *he same in 
essentials but differs in the arrangement of the microwave components. (See the section 
entitled “Experimental Method.”) 

The intrinsic accuracy of the method can be estimated from equations (9) and (10), if 
the calibration errors of the instrument are known. Observational errors are analyzed in 
detail in the section entitled “Experimental Errors.” The maximum errors in e' and e" are 


= ± 


5e n = ± 



!/>■» W I 


2(ja<fa| + |0,rf/3 e |) 


2(|ad0 e |+|/3 e 5a|) 


(13) 


(14) 



where 


1 

6A0 


6AI 

= 


4* 

A <t> — 

A/ 

i 

A/ 2 


5A0 * |S0j| + |60 2 | 
A0 * 0j « 0 2 


( 16 ) 


da = 


5Z 


8A_ 

HAL 

0 


i 


(17) 


The calibration accuracy (maximum values quoted by the manufacturer) of the various com- 
ponents is 

60 (phase shifter) = ± 1/40 rad 
50 (attenuator) » ± 1/60 rad 
8A (phase shifter) * ± 0.1 5 dB (measured value) 

8A attenuator) 0.2 dB 
5AI ~ 2 X 10' s m 


The following a - : representative values for an “as found” soil sample: 

A * 15 dB 
Z. = 12 X 10" 2 m 
A0 ~ 2 rad 

» 

a 15 Np mf 1 
, p e 10 3 rad nr 1 
A/ * 3 X 1 Q“ 3 m 
e' ~ 2.5 X 10“ n Fm" 1 
e" ~ 5.0 X 10~ 13 F nr 1 
v * 37 GHz 

On the basis of these values, 5a = ±0.35 dB and 8(3 c = ±22 rad m~ 1 . For the maximum in- 
trinsic errors in e' and e", maximum errors are 5e' = ±6.5 X 1 0~ 1 3 , 5e'/e / = ±3.0 percent, 
5e" « ±1.9 X 10“ 14 , and 8e"/e" ~ ±3.5 percent. Calculations on a sample of 10 percent 
water content indicate a maximum intrinsic accuracy of be'/e* * ±2.5 percent and 
ge'Ve" = ±2.0 percent. 
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Very lossy samples (water content greater than 10 percent and attenuation greater 
than 30 dB enr 1 ) tend to show considerably larger errors, The primary reason for that is 
the rather short sample length L used in the measurements. Thus a 13.5-pereent sample 
gives Se'/e' ^ *3.5 percent max and 8c n jc” = 14,0 percent max. The same type of sample 
with 20 percent water content shows a maximum error of ±6.0 percent for <=' and ±7,0 per- 
cent for e". A summary of these results and pertinent data are given in table 2, which 
appears later in this document. 


Equations (9) and (10) for the determination of <?' and e" arc correct if the correct 
\ alues of a and p c have been determined. But equation (11) for the determination of a in 
fact determines the amplitude of the wave transmitted by the sample, and neglects reflec- 
tions at the sample interfaces. Further, equation (12) for the determination of the propaga- 
tion constant fi e in effect determines the phase shift of two transmitted waves corresponding 
to two different lengths of sample. It also is an approximation, but a very close one. There- 
fore one must investigate the effect of reflections on the calculation of a and ascertain the 
conditions for proper determination of fi c . 


The problem of finding the transmitted wave for an electromagnetic wave incident on 
a dielectric plate has been solved many times and in many different ways (refs. 13 to 15 ). 
The following discussion briefly sketches the solution of the problem and gives the result, 
using a notation and definitions suitable for the problem. The method is similar to the 
derivation of the Fabry-Perot interferometer equations (ref. 13), with concepts borrowed 
from the solution of Maxwell’s equations. The angle of incidence of the wave is 0° (normal 
incidence), the transverse components of the wave are In the x-y plane, and the direction of 
propagation is in the direction, Let the incident electric wave be A 0x <r 7Z , 7 = a + //3, and 
assume that at z ® 0, fi 0x 6~' tS « 1 . Let R be the complex amplitude reflection coefficient of 
the wave in medium 1 incident on medium 2. R' is the complex amplitude reflection 
coefficient of a wave in medium 2 incident on medium 1 , (See fig. 1 .) The amplitudes of 
the various transmitted and reflected waves can be evaluated at each point from the b jund- 
ary conditions-the continuity of the tangential component of the electric field. For the 
total transmitted amplitude at z~ +L 


£■(« = (1 +/?)(! + R') e - ?£(1 + /?'V 2 t L + R'4 e -*v L + ■■■) (18) 


Summing the geometric series for an infinite number of terms results in 

(1 - R 2 )e~ yL 

£(/)= (19) 

1 - R' 2 e~ 2yL 

If medium 3 is the same as medium 2, R = -R' and 

(1 - R 2 )e~ yL 

£(t) (20) 

1 -J? 2 e’W 
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Figure 1. -Amplitudes of waves transmitted and reflected at sample interfaces. 


For the condition |e~ 27/ 'l /, which can be arranged by making the sample cell length L 

long enough, one can neglect the interference term (1 - R 2 e~ 2yL ), and therefore 

£(') * (1 - R 2 )e^ L 

= (1 ~ Rl) e -(«+m ( 21 ) 

The solution of this equation for ce and the limitation that equation (20) places on the deter- 
mination of p e by equation (12) are discussed in the section entitled “Reflection Errors.” 

EXPERIMENTAL METHOD 


Description of the Microwave Bridge 

Figure 2 shows a schematic diagram of the 37-GHz twin-arm microwave bridge suggested 
by T. Wilheit (personal communication). Attenuator and phase shifter are of the rotary 
vane type. 2 The microwave cw signal was 100-percent square-wave modulated at 1000 Hz 
and detected by a 1N53 crystal and a tuned SWR meter. The microwave source is a back- 
ward wave oscillator with a maximum output of approximately 20 to 25 mW. The un- 
attenuated power reaching the crystal through the reference arm is about 1 1.5 dB down 
from 25 mW, considering losses in isolator's, components, and a 3-dB power division in the 
hybrid T’s. 

The detected noise level output is approximately 0.06 juW at a 1-juW signal level, so 
that the maximum sample attenuation should not exceed 40 dB for a reasonably accurate 
null reading. The noise level is independent of the modulation frequency, indicating that 
the observed noise is probably mostly FM noise rather than 1// crystal noise. Use of a lock- 
in amplifier did not substantially improve the detection sensitivity. 

2 Speciflc equipment was a signal generator 650/659, Alfred Electronics, Inc.; cavity DBD-720-1, Demornay Micro- 
wave Products, Systron-Donner Corp.; precision phase shifter TRG K-528, TRG Division, Control Data Corp.; and a 
precision attenuator HP R382A and SWR meter HP 41 5B, Hewlett-Packard Co. 
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PRECISION 

PHASE 


PRECISION 

ATTENUATOR 


23-dB ISOLATOR, 



There are two hybrid T’s in the bridge, the first used as a power divider, the second as 
a signal comparator. Because proper use of the T’s is of great importance in the operation 
of the bridge, a detailed analysis of their operation is given. 

The properties of the hybrid T are best examined by means of its scattering matrix S 
(ref. 16), The scattering matrix of a T is 

S = (S if ) iff = 1, 2, 3, 4 . (22) 

where S {{ are the scattering coefficients and where subscripts 1 and 2 refer to the two col- 
linear H arms, 3 to the H arm (coplanar with arms 1 and 2), and 4 to the E arm. For match- 
terminated ports, the coefficients have the following properties: s tr s n ’^13 ~* *^23 9 
S 14 = ~»S' 24 , S X1 = $ 22 , and S 34 =0. From S 13 = S 23 it follows then that power into port 3 
divides equally between arms 1 and 2. However, if the reflection coefficients of the loads 
terminating ports 1 and 2 change, the power division will also change. By putting isolators 
in arms 1 and 2 (see fig. 2) the T essentially sees a match-terminated port, and power 
division remains unaffected by the changing of reflection coefficients in the sample arm 


* 
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(due to different samples) and the rather small variations in the reference arm. 3 'File 
important criterion is not equal power division but constancy of power division. The 
power division is a constant at a given frequency but may change with frequency. The 
second hybrid T compares the signals from the reference and sample arm, and the analy- 
sis of its performance is somewhat trickier. There are obviously two conditions that are 
desirable for this T: (1) maximum signal power from both arms 1 and 2 into arm 3, and 
(2) constant power division. It can easily be shown (ref. 16) for a hybrid T that the 
voltage at the detectoi is 


*os-Wi-*n> 


(23 j 


where E n and li n are the input voltages to arms 1 and 2, respectively. This equation is 
true whethu or not the T is matched (i.e., a magic T). If the T is not matched (i.e., a 
hybrid T), then as in the case of the first T, interaction between the T and the load 
(primarily the soil sample) will take place. This interaction will affect the voltage at the 
detector and cause a false reading of magnitude and phase of the wave transmitted by the 
sample. The solution, as in the case of the first T, is to use isolators in arms i and 2. 
One can calculate the power into arm 3 from (ref. 16) 


M 2 + M 2+ M 2 + M 2 = i 


(24) 


Then 


l^ial 2 


S«, 




u\ 


(25) 


i.e., the power into arm 3, which is proportional to \S ] 3 \ 2 may, depending on the magnitude 
ofS l4 , S, 2 , and be considerably less than 1 . In a perfectly matched T, S 12 = S l } ~ 0 
and = 1/2. 

Hence the optimum signal power into arm 3 is 1/2 of the incident power, and conse- 
quently for an unmatched T will be much less. 3 But T’s are never perfectly matched -those 
commercially available have a VSWR in arms 1 and 2 of approximately 1.4 to 1.7; and 
furthermore, the VSWR is fairly frequency dependent. In other words, even in a matched T, 
some interaction between the T and the loads in arms 1 and 2 will take place. A detailed 
theoretical investigation has been made of the effect of load VSWR (or reflection coefficient), 
matched and unmatched T’s, and introduction of isolators, coupled with an experimental 
determination of \S i2 \ , |*S 13 | , \S U \ , and |S n | . Results of this study will be published 
elsewhere. 3 Suffice it to say, the problem was solved by using an unmatched T (see fig. 2), 
and introducing conventional matching devices in arms 3 and 4 to optimize power into arm 
3, or to make |S n | =|5 ?2 | *\S l2 \ =|S 21 | =0. 

Clearly, the matching devices must not destroy the symmetry of the T. (Otherwise the 
T becomes just an asymmetrical four-port device, and the properties of S of equation (22) be- 
come completely irrelevant.) It is easy enough to make an adjustable double stub tuner or 


3 Rcport in preparation. 
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use a commercial stub tuner in H arm 3 to meet the symmetry requirements. The E arm 
tuning element, in order to preserve the symmetry of the T, would have to be introduced 
into the narrow wall of the E arm; i.e., an adjustable inductive post or adjustable H-plane 
stub. It was found, however, that an E-H tuner in the E arm, although asymmetric with 
respect to the symmetry plane of the T, is adequate as a matching element. 

Values in the list following were obtained for the 37-GHz bridge, which shows that 
adequate isolation between 1 and 2 and 2 and 1 was achieved, as well as small and almost 
equal reflections trom the collinear arms; i.e., 1^ J 2 = \$ 22 \ 2 . (All figures are relative to 


Detector arm 3 | 

l s n| 2 

0.01 

fel 2 

0.016 

Reference arm 1 

l«ii| 2 

0.002 

Sample arm 2 

l S 22| 2 

0.0025 

Collinear arms - 

1^121 2 
1 I 

0.00016 


t \ S 2l\ 

0.0001 


Although no measurements were made of the power into arms 1 and 2, a conservative 
estimate can be made of the VSWR corresponding to the relative values of |S 13 | 2 and |5' 11 | 2 , 
The value of the VSWR from arms 1 and 2 is found to be approximately 1.6, comparable 
to commercially available magic T*s. 

The possibility of using hybrid T’s as signal comparators in broadband applications and 
their limitations will be discussed elsewhere. 4 


Sample Cell 

Sample cells were constructed of standard waveguide components; choke-to-choke 
adaptors were used to clamp mica windows to the cover flanges of the test cell waveguide; 
cover-to-cover adaptors provided precision components for changing the length of the cell 
for the measurement of (3 e - A0/A /. Sample cell lengths were held within ±0.0025 cm 
(+0.001 in.), To cover the range of a as a function of moisture content, five sample cells 
were made of lengths 1 1,295, 3.1700, 2.0320, 1.021, and 0.51 1 cm. Figure 3 shows a dia- 


Figure 3,-Microwave specimen cell 
for 37-GHz measurements. 


^Report in preparation. 
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gram of a typical cell. The overall length of the cell was always the same; i.e M 13,7579 cm 
(5.4165 in.) for 37 GHz. held to close tolerances of ±0.0025 cm (^0.001 in.). This allowed 
for easy interchange of all cells in the bridge. 

Sample Preparation and Experimental Procedure 

“As found’’ soil samples held some moisture, depending on ambient humidity. Samples 
were therefore dried. Oven drying, even prolonged heating (20 hr at 393 K (120° €)) 
removed only 2/3 of the moisture (by weight). Vacuum drying was therefore used on all 
samples before any measurements were made or any distilled water was added. Samples were 
vacuum dried for several hours with a small mechanical pump in series with a flask contain- 
ing 500 g of Dryrite. 5 The method of filling the waveguide sample cell was empirical -it 
consisted primarily of vigorous tapping and shaking of the cell as the sample was gradually 
added. Measurement points recorded on the graphs (figs. 6 to 8, shown. later) represent in- 
dependent determinations; i.e., the soil sample was removed from the cell and replaced by a 
fresh sample from the same prepared batch. Relatively large pieces of rock, wood, or clumps 
of soil were eliminated or broken up. AU moisture percentages were by weight based on the 
dry weight of the soil. Moisture measurements were made to within 0.2 to 0.5 percent. 
Maintaining the moisture content while filling the guide is an obvious problem. Extreme 
care was used to expose the moist soil to the air as little as possible. The sample guide was 
sealed off with plastic tape between fillings and during tapping and shaking of the cell. 

Phase angle changes produced by removing the sample from the 3- or 6-mm shims of the 
cell (see fig. 3) were less than 360° to avoid ambiguities in the phase angle readings. At high 
moisture contents the usual procedure was to determine accurately the phase shift for 3 or 
6 mm of sample, and then check this figure against a 1-mm change. A phase change of 360° 
in 1 mm requires a dielectric constant e ' of 64 at 37 GHz-greater than that of pure water. 

To make a direct assessment of the bridge’s accuracy aside from the intrinsic accuracy 
calculations, the specimen cell was replaced by a combination prec : rion attenuator and sliding 
short. This combination permitted simulation of absorption and phase shift. The sample 
arm precision attenuator was changed over a range from 0 to 20 dB at 1-dB intervals. Its 
readings checked the reference arm attenuator to within ±0.1 dB, except two points of 
deviation to ±0.2 dB. The phase errors determined by the micrometer sliding short were ±4°. 
However, the standing wave produced by the 30-dB reflectivity of a 10-dB coupler used in 
this measurement accounts for ±2°. This is in fair agreement with the calibration accuracy 
of ±1.5° of the reference arm rotary phase shifter. 

REFLECTION ERRORS 

An expression for the electromagnetic wave transmission of the sample in the guide was 
derived in the previous section. The measurement of consists in finding the phase shift 


S A registered trademark of the W. A. Hammond Dryrite Co. 
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between two transmitted waves. Let the first wave traverse a length and the second L 2 ; 
then from equation (20) 


(1 -R 2 )e~ yl i 

E (t) - 

1-7?V 2 yl n 


(20a) 



(1 -R 2 )e- yl 2 


1 - K 2 e~ 2yI '2 


(20b) 


These equations show very clearly that the phase change A <j> between Ep and is given by 

A 4>*(L X -L 2 W'-p Q ) (26) 

only if 1 - R 2 e~ 2yL remains constant, which is true if 


A MZL 


Variations in 1 - R 2 e~ 2yL are also minimized by the condition 


e~ 2yL \ «:i 


(27a) 


(27b) 


which must be met to allow calculation of a from equation (21). When A/ ~ L, as in the 
case of high-loss samples, equation (26) is not quite applicable, and 0 e is calculated by an 
iteration procedure together with a from equation (21). 


It remains to find the erreci or 


lUV WIJVVUWIl 


factor” I ~ R 2 of equation (21) on the 


value of a. Before proceeding to a direct solution of equation (21), it will be instructive and 
useful to solve the problem for a given reflection coefficient and calculate the Aa due to 
reflections. It will also be useful to calculate the difference between fi e L and the phase 
<E> = r\L of the transmitted wave. From equation (21) 


(1 ~R 2 )e~(«+M L (28) 

where A is as defined in equation (11). Set A = <r “ >L , and let Aa « a' - a and Aj3 = r? - j3, 
where Aot and A/3 are not necessarily small quantities, and the result is 

! - i? 2 *= e“(A“ + MP)£ (29) 

Writing i? 2 - \R\ 2 e 2 ^° and solving equation (29) for A(3 L and AaL results in 


1 - |/?| 2 cos 2d 

cot Aj 3L = (30) 

\R\ 2 sin 2d 


AaL - - ~ In (1-2 |i?| 2 cos 20) (31) 

Let 1 - \R\ 2 e 2 f° - ce ^ x ; then consideration shows that for 7r/2 < d < 7r and 0 < d < tt/2, 

X varies from 0 <x <tt/ 2 to -*7 t/2<x< 0, respectively. Because x = ~A/3L, then | A/3Z/I <tt/2 
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and IAj 3/J is bounded. Figures 4 and 5 show plots of Aj 3L and Aa/„ respectively, as functions 
of both 1/^i 2 and 0 . They show that even fora relatively large reflection coefficient (I/d 2 ^ 0.3, 
or 1 0 dB),|A/3/,| '■ v 0.3; for representative values of p and /,, A/3//3 1 , or p. This property 

is useful in the solution of equation (28)„ The case for A aL is different AaL increases with- 
out limit* as it obviously must as I/d 2 approaches 1. From the curves, it is seen that Aa is 
comparable or larger than the 8a due to intrinsic errors. Reflection corrections are therefore 
required. 


The solution of equation (28) presents some difficulties. As will be seen in equations 
(32) and (33), and has been seen in equations (9) and ( 10), R can be expressed explicitly as 
R(a, p). Equation (28) is then a complex transcendental equation of degree six in a and p. 

A closed solution of this equation is impossible and it must be solved by iteration. However, 
an approximate calculation of A from a' and p c and substitution in equation (28) with 
<I> ~ p c L will give a very close approximation to a. 

The solution of equation (28) by iteration then should present no difficulties. An 
analytic proof of the convergence of the iteration is presented in appendix A. However, it 
has been assumed so far that all quantities except a are known in equation (28). But d> is 
not. It is difficult to measure without ambiguity because it involves an angle on the order of 
a thousand degrees, but our preceding analysis shows that 4> ~ and is not very sensitive 
to changes in R, (See fig. 5.) It will therefore be possible to obtain an iterated solution of a 
while calculating the correct value of <I> by the same iteration. 


The iteration proceeds as follows: Calculate e' and e" using equations (9) and (10) and 
the approximate values of a and a 1 . A first approximation to R and/?* 1 * is obtained from 
the equation 

Z a V, P c )-Z 0 

/?(D = ■- (32) 

Z (1) («',0„> + 2 n 


where fS c ) is the wave impedance for the TE mode 




(27r/X 0 )(M < A> 1/2 M 1+ /( < */M 

a 2 + p 2 


(33) 


and Z 0 is the characteristic impedance of the guide 



(34) 


All symbols are as defined previously. (The subject of waveguide reflection coefficients is 
treated very obscurely in the literature; only Schelkunoff (ref. 17) gives a fairly clear discus- 
sion, although Goubau (ref. 1 5) is clear once one gets past the notation.) In the calculation 
of R , the wave impedance of the sample was used rather than its input load impedance Z L , 
a justifiable procedure because there is always enough attenuation in the sample so that 
Z(a, P) approximates Z L very closely. The first approximation to R is now substituted in 
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0, PHASE ANGLE OF REFLECTION COEFFICIENT, deg 


Figure 5,~Attenuation constant versus pfiase angle of reflection coefficient 
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equation 1 28) for a and <I> to obtain a second approximation to a, <* ,2) and a first approxima- 
tion to <Ik «I» M \ The general iteration equation is clearly 

( 1 - \R\M2 c 2jO Ui )( r^* = r «'I, ( ,jJ ( L ( 35 ) 


The expression on the right remains constant throughout the iteration. 


For relatively low-lor samples (moisture content less than about 10 percent), direct 
determine ';f 0 t was possible. The conditions for this determination were given in equa- 
tions (2 7a) d't: (27b). iiigh-loss samples require a different procedure, because it is no 
longer possible to satisfy equation (27a). In fact, it becomes more convenient and straight- 
forward simply to determine the phase shift of the transmitted wave in kmgth L of the sam- 
ple aim then determine j3 f from equation (28). Two measurements are made as ueforc, the 
first measurement w Ji the sample in the cell of length L. The equation for the transmitted 
wave is then, from equation (28), 

( 1 - R 2 )<r (a+ ^ L = cr *’ 1 (T hl (36) 


A second measurement without the sample gives e~ // 3 o L for the transmitted beam. If the 
phase determination for the first measurement gives <1^ , <I > 2 for the second, and A<F = 'J>. - <I> 0 , 
then dearly, from equation (36) 


AT> = rjl ~ P Q L 
A <p 27 r 

1 \‘A 

In other words, in equation (36) the quantities a! and r\ are now known and it remains to 
find a and p g . From the preceding analysis, it was found that p ~ fi e . In fact, from figure 4 
one can estimate for/? ~-0,5 /170° that A/3 = 50 rad m ~ 1 if L = 3 X 10 " 3 m. One can 
therefore formulate the iteration equation for a and p e by writing equation (35) in the form 

(1 - 0») e ~ u<i+m e-^ i+l)L 

and setting |3 (1) ~ 77 and 0 ^ = a r . 


The condition for convergence of the iteration equation is given by 

2RdR 
L da 


(37) 


Equation (37) is derived in appendix A. For a very lossy sample (15 percent water content) 


and 


dR 

00! 


- 10- 5 


L> 


2 R 


9 R 
da 


= 1 X 10 - 5 m 


Hence convergence is assured for even very short sample lengths. 
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Table 1 is a reproduction of a computer printout of successive iterations for samples 
witli various water contents. As expected, the first iteration results in values of a, e', and e n 
(and for very lossy samples) that differ little from subsequent iterations. R' and R" are 
the real and imaginary parts of the reflection coefficient R ~ R' + jR", 


Reflection corrections for relatively low-loss samples (less than 10 percent water) are 
generally small compared to intrinsic and observational errors. The real part of e is almost 
unaffected by a t 0 being very much larger than a. Hence corrections to e' are too small to be 
of importance in this low water content range. For the imaginary part of e, which is directly 
proportional to a, the correction can amount to several percent. 


Very lossy samples (moisture content over 10 percent) require correction in a as well as 
as shown previously. The constants a and are now of the same order of magnitude, 
and, consequently, there will l v ; substantial corrections of both e ; and e". As expected, the 
measured (x (i,e„ a') is much too high, and the measured dielectric constant e', too low. 

Oddly enough, the correction of e" is less than expected because the decrease in a is partially 
compensated by an increase in . Reflection corrections for a and are approximately 10 
and 1 percent, respectively. 


EXPERIMENTAL ERRORS 

The intrinsic accuracy of the method, discussed previously, k the maximum accuracy 
obtainable considering calibration errors of attenuators, phase shifters, and measurement 
errors in the cell lengths (and shim lengths). In this section, the error calculations are con- 
cluded by considering in addition the observational errors. (Errors in the determination of 
A 0 and A^ are not considered because they are very small compared to all other errors.) 
This classification of errors is somewhat arbitrary but useful. Observational errors in the 
reading of A -and 0 change with the changing ratio of detected signal level to detected noise 
in the crystal null detector. Intrinsic and observational errors are considered statistically 
independent and are combined by calculating the rms error in e from 



where 6 2 a ~ 6 2 <*j + 8 2 ot 2 and where dot^ is the estimated intrinsic and 8a 2 the estimated 
mean observational error. The process is similar for 0. 

The calculation of 5e rms is rather tedious and its details are omitted. Calculations were 
made for samples with representative moisture contents (i.e., 0.6, 10, 13.5, and 20 percent). 
Errors for other moisture concentrations were assumed to be between the appropriate limits. 
Table 2 summarizes the results and shows both maximum intrinsic errors and total rms 
errors. Graphs of e versus percent of water show scatter that in many cases far exceeds the 
rms error. This experimental point scatter v% clearly caused by difficulties of getting repro- 
ducible samples and reproducible sample pa Jang. 
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Table 2. Representative Errors and Correction Factors 


Sample 

Maximum 
intrinsic 
error, a 
percent 

Total 
ms error, 
percent 

Reflection j 
correction 
percent 


be 1 

be" 

be img 

be" 

oc rmi 

be' 

be" 

M5, as found, 0.66 percent ICO, 
a ® 15 Np nr 1 , ft 1 Q 3 rad in" I , 

L ~ 12 X 10 2 m, A0 ~ 2 rad, 
A/~3X 10 ’ 3 m, c' = 2.4 X 10 11 , 
e" = 5.5 X 10 13 

±3.0 

+2.5 

±3.6 

±3.5 

0,0 

-6 

M5, 10 percent 300 Np nfl, 

1500 rad in" 1 ,// ^T. 5 X 10~ 2 m, 
A/ ~ 3 X 10" 3 m, A0 ~ 3 rad, 

(f “ 4.0 X 10" n ,c" = 1.5 X1Q~ U 

±2.0 

±2.0 

±3.5 

±2.5 

“0.3 

-6 

F2, 13.5 percent H 2 0,a~ 440 Np nr 1 , 
/3 ~ 1900 rad nr 1 , L ~ 5 X 1 0 ' 3 m, 
A0 ~ 6.5 rad, e ' = 5.5 X 10” 1 1 , 
e" = 2.5 X 10” 11 

±3.3 

±4.0 

±3.3 

±3.0 

+3.5 

“12 

F2, 20 percent H 2 O, a ~ 700 Np m -1 , 

(3 ~ 2300 radm”l,£ ~3 X 10" 3 m, 

Atj> ~ 5.3 rad, c' = 7,3 X 10" 11 , 

= 10-11 
» Ik/ f 1 JL V 

±6.0 

±7.0 

±5.6 

±6.0 

+5 

-10 


fl Scc '•.cction entitled “Method of Determining a and /? c .” 
b See section entitled “Reflection Errors,” 


DISCUSSION OF RESULTS 

Published data on dielectric constants of soils with varying water contents are somewhat 
meager. The first published results of soils with a number of different moisture percentages 
in the rf range (at 10 MHz) seem to have been those of Smith-Rose (ref. 18). Some of the 
first microwave measurements apparently were made by Ford and Oliver (ref, 19) at 9 cm on 
dry and “very wet” soil. They used a combination of in situ reflection and absorption be- 
tween receiver and generator antennas. More recent are the X-band measurements on dry 
soils, and soils with 6, 10, and 25 percent water by Straiton and Tolbert (ref. 5) and by von 
Hippel (ref. 20) on 2, 4, and 16 percent water-content samples. There is a solitary measure- 
ment at 19 GHz on dry sand by Hertel, Straiton, and Tolbert (ref. 10). 

Interestingly enough, dry soil samples taken from diverse and obscure locations and of 
various compositions show a surprising degree of agreement in the value of e. Measurements 
by Straiton and Tolbert (ref. 5) show a variation from 2.5 to 2.85 X 10" n Fm" 1 , von 
Hippel’s (ref. 20) from 1.9 to 2.2 X 10" 11 , and this study at 37 GHz for three soil samples 
gave values from 2,5 to 2.7 X 10" 11 . However, the recent results by Wiebe (ref. 21) at 
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X-band on soil samples from Texas show dry soil values substantially higher than the above, 
some as high as 3.9 X 10" 1 1 . 

With added water, there is a radical departure of values found in this study from those 
in published data, Dielectric constants e' ranged from 3.2 to 4.5 X 10" 11 F nr 1 at 10 per- 
cent water; the corresponding values of Wiebe show a spread from 4.5 to 8 X 10" 1 1 . At a 
moisture content of 30 percent, the disparity becomes even greater. The values found in this 
study rise to approximately 9 X 10” 11 ; however, Wiebe’s data show a rise to 28 X 10" 11 , 
with a minimum of roughly 16 X 10" u . The data by Straiton and Tolbert (ref. 5) show a 
similar strong increase in e f as water is added. 

A partial explanation of this behavior can be attributed to the variation of the dielectric 
const m t of water with frequency (ref. 12, fig. 1 1.7). The real part of e decreases by a factor 
of when the frequency is changed from 9 to 37 GHz; the imaginary part e", on the other 
hand, changes relatively little in the same range. The latter behavior seems to be reflected in 
the data for e". Values found in this study for 30 percent moisture content are in the 
neighborhood of 6 X 10" 11 ; Wiebe’s values are approximately 8 X 10" 11 . 

Results of this study for e' and e" for soil samples of loamy find sand (sample designa- 
tion M5), fine sandy loam (L3), and sandy clay loam (F2) are plotted as functions of percent 
water by dry soil weight in figures 6 to 8. 



figure 6.— Values of e sind e” versus 
water content. Loamy fine sand sample 
M5, 37 GHz. 


20 




Figure 7 .-Values of e' and e" versus water content 
Sandy clay loam sample F2, 37 GHz. 


Figure 8.-Values of e' and e” versus water content. 
Fine r>andy loam sample L3, 37 GHz. 



0 5 10 16 ,'Q 

HjO (percent by dry wt) 
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Table 3, -Textural Analysis of Phoenix, Ariz., Area Soils 


Sample designation 

Soil texture 

Percentage composition, dry 

Sand 

Silt 

Clay 

M5 

Loamy fine sand 

88.0 

7.3 

4.7 

L3 

Fine sandy loam 

48.0 

34.0 

18.0 

F2 

Sandy day loam 

56.0 

26.7 

17.3 


Origin of the samples and their textural analysis are given in table 3. A!1 dielectric con- 
stant determinations are shown independently on the graphs and no averages are taken. The 
reason for this is simply to exhibit the scatter of the results caused by differences in sample 
packing and sample variability. Measurements at approximately 0.6 percent water correspond 
to “as found” soil samples. 

The most interesting and unexpected feature of the curves of e 1 versus percent water is 
the behavior of e' at relatively low moisture content (up to approximately 10 percent). Con- 
sidering the extremely high dielectric constant of water, an immediate and substantial in- 
crease in the dielectric constant of the mixture was expected as water was added to the soil. 
This expectation was also based on the published data (refs. 5 and 21). But the dielectric 
constant observed was almost independent of or decreasing with moisture content, reaching 
a broad minimum in the neighborhood of 6 percent. 

Further increase in water content resulted in a steep increase of e'. The minimum seems 
to be preceded by a broad maximum-however, the data are not good enough to demonstrate 
this clearly. A similar drop in e' is evident in data published by von Hippel (ref. 20) and 
Bdgerton (ref. 22). However, Wiebe’s data (ref. 21) show a smooth increase of e' with water 
content; that trend may have been the result of the sample packing method (a controlled 
impact hammer) and the lack of detailed data in the 0- to 10-percent moisture range. 

On the whole, the results for three soil samples show good agreement. The maximum 
dielectric constant e f falls in the range from 8 to 9 X 10“ u F m" 1 and e" from 5.5 to 
7 X 10“ 11 F nr 1 . The initial variation of e' with moisture is relatively small, about ±30 per- 
cent, but above 10 percent water, e' increases at the rate of about 1 X 10“ 11 Fm" 1 per per- 
cent of water. This extremely rapid increase of e' and e" with moisture is not reflected in 
the free-space reflection coefficient; as a result, the microwave brightness temperature 
(eq. (1)) variation with water content still has manageable proportions. In fact, one can 
easily derive an expression for the change in brightness temperature as a function of mois- 
ture. The details of the calculation are giyen in appendix B. 

If the soil sample has a dielectric constant e « e' - /e", then the free-space amplitude 
reflection coefficient for normal incidence is 

1 - deb 

R = (39) 

1 + dd* 
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where 



The power reflection coefficient RR * becomes 



1 - 2d cos v? + d 2 

RR * = 

1 + 2c? cos ip + d 2 

but cos v? ~ t , and 


( 40 ) 


RR> - (i^J (4.) 

Calculating the total differential of RR *, and assuming d(RR*)/dy ~ 0, one ’finds 

d(RR'*) 2 {d 2 - 1) / de' de"\ 

= + U — + e " (42) 

dx egc/ 3 (J+l) 4 \ dx dx) 

where x is the percentage moisture. Because the radiometric brightness temperature is given 
by T b = (1 - RR*)T S0 £, the result is 

dT h d(RR*) 

1 = T (AVi 

, , 1 sofi 

dx dx 


For sandy clay loam (sample F2) with 13 percent water, e'/e e ~ 5, e"/e Q ~ 1.5 (from fig. 
7), and 


d(RR*) 

dx 


de' 

0.03— +0.01 

dx 


de" 

dx 


From equation (43) 

dT b ~ 

- - 12 K/percent moisture 

dx 

for 

de' „ de" „ 4 t . 

— = =1 F m" 1 /percent moisture 

dx dx 


CONCLUSIONS 

Measurements have been made of the complex dielectric constant of three types of 
soil from the Phoenix, Ariz., area at 37 GHz. From 0 to 40 percent distilled water v/as 
added to the dried soils in more or less regular percentage increments. Special care was used 
to establish the relative insensitivity of e' to water content up tc 10 percent and the steep 
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increase in the 10- to 20-percent region. The usefulness of the dielectric data for thermal 
sensing would therefore be limited to approximately the 10- to 20-percent moisture region. 
In that region the microwave brightness temperature change per percent moisture change in 
the soil is approximately - 1 2 K per percent change. 


All measurements were corrected for reflections at the sample interfaces. These cor- 
rections primarily affect the values of a and e " at a moisture content of up to 10 percent. 

For high-loss samples, corrections must be applied to both a and (3. The error in a at 10 per- 
cent water content is approximately 5 percent, and at 30 percent water content, 9 percent 
in a and 3 percent in (3. Graphs were drawn to allow the corrections Aot/ot and A/3//3 to be 
estimated for a given reflection coefficient and length of sample. Actual calculation of a, 

13, e\ and e " was done by machine iteration of transcendental equation (35). Two experi- 
mental methods were evolved to deal with lossy samples. Low-loss sample phase shifts were 
determined from a conveniently small section of the sample cell; when losses exceeded about 
50 dB cm" 1 , phase shifts were measured directly for the total length of the cell. Accordingly, 
two different computer programs were used in the calculation of the dielectric constants. 


Limit expressions were derived for the proper experimental determination of (3 e , 
suppression of interference effects, and the convergence of the iteration calculation. 


the 


Although reflection errors can (see analysis in section entitled “Reflection Errors”) be 
reduced to zero if the sample cell is long enough (i.e., if A odL is constant and Aa 0 as 
l-+a)\ however, this is not a practical approach. The limited microwave power and sensitiv- 
ity of the detection equipment limit the total amount of attenuation that can still be 
detected long before Act and A/3 0. (See discussion in section entitled “Experimental 

Method.”) 
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Appendix A 


PROOF OF CONVERGENCE OF ITERATION 


Equation (28) can be rewritten in the form 

[ 1 - R 2 (a, 0)]e- < "V 4 ' = c® 1 e~ aL (A-l) 

For a given measurement, the expression on the right of the equal sign is a constant. In the 
iteration calculation both a and <E> are iterated, but the change of $ with R is small, as dis- 
cussed in the section entitled “Reflection Errors.” Assume for the sake of simplicity that 
is constant. Then equation (A-l) can be rewritten in the general functional form 

f(a)F(oi) = F(cl q ) (A-2) 

where 

f(a) = 1 - R 2 (a) 

F(a) = e~ aL e^ 

and « 0 is the experimentally obtained data and represents the first approximation to a. The 
term F(ot Q ) is a constant. 

Write the iteration equation in the form 

f(« r )F(ot r+l ) = F(a Q ) (A-3) 

and let 


a = a r + 8a r 


Expanding both/(a r ) and F(oi r+1 ) by Taylor’s theorem and using f(oi)F(ct) = F(ce Q ) gives 

9 f/da F(a) 


&Ci — 

r+1 dF/da f(ci) 


(A-4) 


Let a 0 - aj = 8o: 0 , then 


8a 


; r+1 


8 f/da F(a) ' 
8F/8a /(oi) 


\\ 

) K 


(A-5) 


or 


^ a r + 1 0 
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From equation (A-2), it is seen that 


One can write 


Then 


and from 


/’ 


2R 9 R 

F' f 


L 9a 


R=R[2; 2 (a)] 


m 

9/^ 9Z 2 

9 a 

3 Z 2 3a 


2Z 0 

9Z 

(Z 2 +Z 0 )' 


M o w0 

Z. 

3 

2 

a 2 + P 2 


(the imaginary part can be neglected) 

dZ 2 2/jIq coce 
3a (3 3 


With 

a ~ 200 Np m~ 1 
P ~ 2000 rad m" 1 
UM 0 ~29X 10 4 
Z 0 ~ 460 n 
z 2 ~ 200 n 


3 R 4j a 0 coaZ 0 
3a (Z 2 +Z 0 ) 2 j3 


the result is 



Appendix B 

CHANGE OF BRIGHTNESS TEMPERATURE WITH MOISTURE 


The propagation constant for a transverse electromagnetic (TEM) wave (ref, 12, eq. 
(16)) is given by 


Make the substitution 


and let o - 0, then 


7 = VJwfr + /coe) 


_ _ j ,>n 

e = e - je 


(B-l) 


7 = y(/coju(cue'' +/<oe / ) 


(B-2) 


Comparison with equation (B-l) shows that one can make the formal substitution a = cue". 


mi /*. 


i iie free wave impedance uf a medium is defined (ref, 12, eq, (17)) by 


Z = 


/cu/x 


o + /cue 

From the above it follows that for a lossy medium of e « e' - /e”, 


• (B-3) 


Z = 


7M 


e" + /e' 


(B-4) 


The amplitude reflection coefficient for normal incidence for the E vector of a TEM wave 
incident on a medium of wave impedance Z 2 (ref. 17) is 

Z« ~ Z.\ 


r e 


Z 2 +Z 0 


(B-5) 


Substitution for Z 2 and Z 0 = gives 




1 - s/e'/e 0 - /e"/e c 
1 + \/e'le 0 - /e"/e c 


(B-6) 
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Let 


then 


And one can write 



V t o 


e'/e„ 

cos IV = " — 

\/V/e 0 ) 2 *He"/e 0 ) 2 



4 



1/4 


1 - deb 

R e = T (B-7) 

1 + de’* 

The power reflection coefficient then becomes 

****■(&$ ,m > 

Because e' > e", cos 2<p — 1 , andcos^~ 1. The total differential o(RR* is calculated from the 
partial differentials, bR/bd, 3 d/de', and bd/be", resulting in 

2 (d 2 - l) 

d(R E R%) = (e' de 1 + e " de ") ( B-9) 

6q(1 + d)*d 3 
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